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We generalize the mean-field theory for the spinless Bose-Hubbard model to account for the 
different types of superfluid phases that can arise in the spin-1 case. In particular, our mean-field 
theory can distinguish polar and ferromagnetic superfluids, Mott insulators which arise at integer 
fillings at zero temperature, and normal Bose liquids into which the Mott insulators evolve at finite 
temperatures. We find, in contrast to the spinless case, that several of the superfluid-Mott insulator 
transitions are first-order at finite temperatures. Our systematic study yields rich phase diagrams 
that include, first-order and second-order transitions, and a variety of tricritical points. We discuss 
the possibility of realizing such phase diagrams in experimental systems. 
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I. INTRODUCTION 

Experimental investigations of ultracold atoms in op- 
tical lattices have opened up a new realm in the study 
of quantum phase transitions (QPT) 1 - 2 . The super- 
fluid (SF) to Mott-insulator (MI) transition has been ob- 
served in spin-polarized 87 Rb atoms trapped in a three- 
dimensional, optical-lattice potential 3 -, by changing the 
strength of the onsite potential, as predicted theoretically 
by studies of the spinless Bose-Hubbard model-' 5 . Fur- 
thermore, technical advances in the trapping of atoms by 
purely optical means* have enhanced the interest in the 
study of quantum magnetism in confined dilute atomic 
gases. Alkali atoms with nuclear spin / = 3/2, such as 
23 N a, 39 K, and 87 Rb, have hyperfine spin F = 1. In 
conventional magnetic traps, these spins are frozen, so 
the atoms can be treated as spinless bosons; by contrast, 
in purely optical traps, these spins are free, so the Bose 
condensates, which form at low temperatures, can have a 
spinor nature 7 -^ and the SF-MI transition can be modi- 
fied 9 - ' 10 : 11 . In the spinless case, the SF-MI transitions are 
controlled by the interaction Uq between bosons at the 
same site. As Uo increases beyond a critical value, the SF 
phase undergoes a continuous transition to an MI phase 
in which the number of bosons at every site is an integer. 
This transition is reflected in the development of a gap at 
the transition. When the spin is nonzero such a gap also 
develops at SF-MI transitions, but the properties of the 
phases and the natures of these transitions are modified 
by the spin degrees of freedom. 

Theoretical work on this problem has dealt primar- 
ily with the properties of spinor condensates by using 
a continuum, effective, low-energy Hamiltonian. Such a 
Hamiltonian suffices if one is interested in the natures of 



the superfluid phases, which can be polar or ferromag- 
netic, and in their excitations, which include vector or 
quadrupolar spin waves and topological defects 7 -^. How- 
ever, if we want to study the SF-MI transitions we must 
use a lattice model such as the spin-one Bose-Hubbard 
model. Some group a 9 ' 10 ' 11 have initiated such an investi- 
gation by obtaining the zero-temperature phase diagram 
of this model in a mean-field approximation. The topol- 
ogy of this phase diagram for the spin-1 Bose-Hubbard 
model is similar to that of its spinless counterpart; but, 
in the spin-1 model, the superfluid phases can be either 
polar or ferromagnetic depending on whether the spin- 
dependent interaction favors or disfavors the formation 
of singlets. In the former case the SF-MI phase transi- 
tion is continuous, if the density of bosons per site p is 
an odd number, but first-order, if p is an even number. 

We have two main goals in this paper: The first is 
to give a global view of the zero-temperature, mean- 
field-theory phase diagram of the spin-1 Bose Hubbard 
model emphasizing issues of first-order coexistence that 
have not been highlighted so far. The second is to 
generalize this mean-field theory to finite temperatures 
T > and thus obtain the finite-temperature, mean- 
field-theory phase diagram for this model. 

The spin-1 Bose-Hubbard model is defined by the 
Hamiltonian 



<i,j>,<r i 

+yE(^ 2 - 2 ^)-E^. (!) 

i i 

where the first is the kinetic energy associated with the 
hopping of bosons between nearest-neighbor pairs of sites 
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< i,j > with amplitude t, a] a {a^ a ) is the boson creation 
(annihilation) operator at site i with spin component a 
(which can assume the values 1,0,-1), rii a = a ia .ai tCT ; 

tii = Y.<j n i<y and Fi = a ' c H,a F T,"' a i,<'' are ' respec- 
tively, the total boson number and spin on site i, and 

F aM i are the standard spin-one matrices; Uo is the on- 
site Hubbard repulsion and U 2 the energy for nonzero 
spin configurations on a site. The origin of such a 
spin-dependent term lies in the difference between the 
scattering lengths do and a 2 , for S = and S = 2 
channels^, respectively; in terms of these lengths Uq = 
4:Trh 2 (a + 2a 2 )/3M and U 2 = ^h 2 (a 2 - a )/3M, where 
M is the mass of the atomi. For 23 Na, a 2 — 54.7as and 
a = 49.4as, where as is the Bohr radius, so U 2 > 0, 
whereas for s7 Rb, a 2 = (107±4)ob and a = (110±4)ob, 
so U 2 can be negative. The parabolic trapping potential 
with strength Vr is represented by the site-dependent 
chemical potential fa = /i ~ Vr\Ri\ 2 , where Ri is the 
distance of site i from the center of the trap and /1 is a 
uniform chemical potential that controls the mean den- 
sity of the bosons. In this study we neglect the trap 
potential (i.e., we set Vr — 0) and focus on the effects of 
the spin degrees of freedom. 

The zero-temperature phase diagram of model (Q} has 
been obtained in the mean-field approximation by some 
group o 9 ' 10 ' 11 . We have extended these studies signifi- 
cantly Before presenting the details of our work, we give 
a qualitative overview of our new results. 

Consider first the case U 2 — 0: We might expect the 
spin-1 and spinless Bose- Hubbard models to have same 
phase diagram in this case since the ground-state en- 
ergy does not depend on the spin. This is superficially 
true at T = in so far as the SF-MI phase bound- 
aries for both these models overlap. However, as we will 
show, the SF phase is highly degenerate in the spin-1 
case; and for T > the SF-MI transition becomes first- 
order and, eventually, continuous again. Thus the finite- 
temperature phase diagram has a rich topology with first- 
order boundaries evolving into continuous ones at tricrit- 
ical points. 

If U 2 the onsite interaction between the bosons 
becomes spin dependent. It turns out that we must dis- 
tinguish between the cases U 2 < and U 2 > 0. The 
former yield a phase diagram that is very similar to the 
one for U 2 — 0; the major qualitative difference arises in 
the nature of the SF phase that is now a ferromagnetic 
superfluid. 

There are many differences between the phase dia- 
grams of the spin-1 model with U 2 = and U 2 > 0. 
If U 2 > the SF phase is a polar superfluid. Further- 
more, even at T = 0, the SF-MI transitions are different 
for odd and even densities. For odd densities, the T = 
SF-MI transition is continuous as for U 2 = 0; however, 
for even densities, this SF-MI transition turns out to be 
first-order because of the formation of singlets that also 
stabilize the MI phase considerably. 

At finite temperatures the MI phases evolve without 



a singularity into a normal Bose liquid (NBL). These 
are really not distinct phases but, as we will show, the 
compressibility k can be used effectively to delineate the 
crossover between MI and NBL regions. 

To present our results in detail we must introduce our 
mean-field theory. We do this in Sec. II. Our results are 
given in Sec. III. We end with a discussion in Sec. IV. 

II. MEAN-FIELD THEORY 

Mean-field theory has been very successful in obtaining 
the phase diagram for the spinless Bose-Hubbard model. 
There are three formulations of this mean-field theory: 
one uses a model with infinite-range interactions, another 
a Gutzwiller-type wave function, and a third 5 , which we 
follow, a decoupling approximation. The unique fea- 
ture of this decoupling scheme is that, unlike conven- 
tional mean-field theories, it does not decouple the in- 
teraction term to obtain an effective, one-particle prob- 
lem but, instead, decouples the hopping term to obtain 
an effective, one-site problem. This one-site problem 
is then solved self-consistently. We generalize this de- 
coupling procedure to the spin-1 case as follows 5 -: In 
the identity a\ a a ha = (a\ a - (o| )<7 ))(a 3 > - (a jt<T )) + 

( a t<r) a 3f+ a \,*( a o>°) ~ ( a L)( a ^)> where (O) denotes 
the equilibrium value of an operator O, we neglect the 
first term that is quadratic in deviations from the equi- 
librium value. Thus 

a \,a a h<y - ( a l*) a j,° + a i,a( a j,<r) 

-(alj(a^); (2) 

since we expect superfluid phases, it is natural to intro- 
duce the superfluid order parameters 

= (4,*) = ( a '^)> ( 3 ) 

for a = 1, 0, —1. We consider equilibrium states with uni- 
form phases, so we choose these order parameters to be 
real. Given the decoupling approximation {2J the Hamil- 
tonian ([1]) can be written as a sum over single-site, mean- 
field Hamiltonians: 

w=5>r, (4) 

i 

where 

Uf F = l ^h i {n i -l)+ l ^{F?-2n i )-»n i 

-Mala+^,a)+Y,\^\ 2 - ( 5 ) 
a 

We set the energy scale by choosing zt = 1, where z is the 
number of nearest neighbors. At least one of the order 
parameters ip a is nonzero in a superfluid phase. In order 
to calculate i/v m our mean-field theory, we first obtain 
the matrix elements of the mean-field Hamiltonian Ti.f IF 
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in the onsite, occupation-number basis {|n_i,no,ni >} 
truncated at a finite value n m ax of the total number of 
bosons per site n = $3 CT ? V- I n most of our studies we 
use n m ax = 4 for which the mean-field Hamiltonian is a 
36 x 36 matrb*i£. We diagonalize this matrix to obtain 
its eigenvalues E a and eigenvectors | ip a }- 



Hf*\<p a )=£ a \<p a ); 



(6) 



we suppress the site index i on these eigenvalues and 
eigenvectors since all the phases we consider are spatially 
uniform. 

We now obtain the variational the free energy 
T(n,U ,U 2 ,T;i; a ) = -TlnZ(fi,U ,U 2 ,T;i; a ), (7) 
where Z(p, Uq, U 2 , T; i/v) is the partition function 



Z(p,U ,U 2 ,T;ip (T ) 



(8) 



with the Boltzmann constant ks chosen to be 1. The 
variational free energy T must be minimized with re- 
spect to the order parameters ip a , i.e., we must solve the 
equations dJ-/dtp a = for a — 1, 0, —1. These equations 
can be recast as self-consistency conditions for so- 
lutions of these self-consistency conditions correspond to 
extrema of T . In case there is more than one solution, 
we must pick the one that yields the global minimum of 
T . [At a first-order phase boundary T has two, equally 
deep, global minima.] The values of ip a and T at the 
global minimum yield the equilibrium order parameters 
and free energy J- eq . In our mean-field theory, the 
superfluid density is 



Lie 



q |2 



(9) 



The magnetic properties of the superfluid phases of this 
model are obtained from — 



(F) = 



ell 



e 9|2 



(10) 



the explicit forms of the spin-1 matrices now yield 

1,2 



(F) 
(F) 2 = 2 



A 0M>o + V>-iV>o) . (^-fe) 



E.h/vl 2 
(EJ^I 2 ) 2 



EJVvl 2 
(V* 2 -^) 2 
(EJVvl 2 ) 2 ' 



(ii) 



where x and z are unit vectors in spin space and we 
suppress the superscript eq for notational convenience; all 
ip a used here and henceforth are actually ip% q . Superfluid 
states with (F) = and (F) 2 = 1 are referred to as polar 
and ferromagnetic, respectively. The polar stated has an 
order-parameter manifold (U(l) x S 2 )/Z 2 , where U(l) 
denotes the phase angle 9, S 2 refers to the directions n 
on the surface of a unit sphere (on which orientations 
are specified by the angles (a, (3) of the spin quantization 



axis), and Z2 arises because of the symmetry of this state 
under the simultaneous transformations 8 — > 6 + 7r and 
n — > — n. Thus the superfluid order parameters can be 
written as 




■^e~ la sin 13' 



lp s e 



cos (3 



(12) 



! sin f3. 



Similarly, since the ferromagnetic superfluid state has 
an order-parameter manifold^ with the symmetry group 
SO(3), 




= VP s e 



e~ la cos 2 f 
\pi cos ^ sin 4 

m „• 2 j9 



(13) 



where a, (3, and r are Euler angles. 

We consider only spatially uniform supcrfluids in equi- 
librium, so it suffices to use real order parameters. Thus, 
for the polar superfluid we have the following possibili- 
ties: (i) -01 = ip-i > and ijj — with 9 = a = = n/2 
or 9 = —a = (3 = ir/2; and (ii) tpi = 4>-i = and t/)Q > 
with (3 = or ir, 8 = or 7r, and < a < 2ir. Similarly 
for the ferromagnetic superfluid tpi = ip—i, [3 — n/2, 
a = 0, < 6 = t < 2tt, and ^ = V^ipi. 

The equilibrium density p and compressibility k can 
be obtained from 

P = -^- = ^T,e-^(^\h\^), (14) 

where E a and | 
mum of T and 



are £ a and | tp a ) at the global mini- 



dp 

dp 



(15) 



The three quantities p s , (F) 2 , and n together deter- 
mine the thermodynamic phase of model (1) for any 
point in the parameter space {p, Uq, U 2 , T} as given in 
the Table HI Strictly speaking there is no distinction be- 
tween the Mott insulator (MI) and the normal Bose liq- 
uid (NBL); the former exists at T = and has k = 0; it 
evolves without any singularity into the NBL at T > 0; 
at low T, the compressibility k is exponentially small in 
the NBL so one can think of it as an MI phase; at high 
T, where k is substantially different from 0, it is best to 
think of this phase as a normal Bose liquid. It is con- 
venient, therefore, to define a crossover boundary above 
which k is substantial; we use the criterion k = kx = 0.02 
to obtain the MI-NBL crossover boundary shown in some 
of our phase diagrams. We must remember of course that 
this is not a strict phase boundary and it depends on the 
value we choose for the crossover compressibility kx ■ 



III. RESULTS 

We are now in a position to present the results of our 
mean-field theory. It is necessary to distinguish between 
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FIG. 1: (a) Phase diagrams in the (fi — Uo) plane for U2 = 0: Solid lines indicate the T = continuous phase boundaries 
between SF and MI phases; these phase boundaries evolve into first-order boundaries at finite (but low) temperatures as shown 
by the representative dashed lines for T = 0.05; at higher temperatures these first-order boundaries become continuous again 
at lines of tricritical points. Pseudo-grayscale plots at T = of the variational ground-state energy £0 for (b) U2 = 0, (c) 
U2/U0 = 0.03, and (d) U2/U0 = —0.03, respectively; the four degenerate minima in (c) and (d) show that the SF phase is polar 
in the former and ferromagnetic in the latter; in case (b), i.e., XJi = 0, the SF phase is infinitely degenerate (see text). 



5 




FIG. 2: Representative plots of (a) p s and (b) p versus p for Uo = 4, 5, 6, U% = at T = 0. Similar plots at T = 0.05 are given 
in (c) and (d). In the MI phases p s = and p is an integer [= 1 in (b) for Uo = 6] for T = 0; for < T the MI phase evolves 
without a singularity (see text) into the normal Bose liquid (NBL) in which p is exponentially close to an integer [= 1 in (d) for 
Uo = 5, 6]. As the temperature increases from zero, the MI phases grow at the expense of SF phase and the SF-MI transition 
becomes first order [see Fig.([T]. 



three qualitatively different regimes: (1) U2/U0 — 0; (2) 
U2/U0 > (we use U2/U0 — 0.03 since this is appropriate 
for 23 No); and (3) U 2 /U < 0, as in 87 Kb (for specificity 
we use U2/U0 = -0.03). 

We first consider U2/U0 — 0, which can be achieved 
when the scattering lengths are equal, i.e., do = 02- In 
this case the onsite interaction between bosons is spin in- 
dependent. This leads to an infinitely degenerate super- 
fluid state: Specifically, for a given value of the superfiuid 
density p s , the three order parameters i\) a ,o — — 1 } 0, 1, 
can have any magnitudes that satisfy Eq. ©; e.g., if 
we make the specific choice ip-i — ipx, then the pseudo- 
grayscale plot of Fig. Q] (b) shows that the minima of the 
variational mean-field energy at T = lie on the ellipse 
2ip1 + -0q = p s . This degeneracy makes the superfiuid 
phase of the spin-1 Bose-Hubbard model different from 



its spin-0 counterpart even if U2 — 0; and it implies that 
an infinite number of SF phases coexist at U2 = 0. How- 
ever, the zero-temperature phase diagram of Fig. Q] (a) is 
the same as that of the spinless Bose-Hubbard model^; 
and, in particular, lobes of the MI phase are separated 
from the SF phase by the SF-MI boundaries that are all 
continuous at T — 0; the density p is fixed at integral 
values in each MI lobe. 

Striking differences between the spin-1 and spinless 
cases appear at finite temperatures. We demonstrate this 
in Figs.[]Ja) and[21 the former compares phase diagrams 
at T — and T = 0.05 and the latter presents plots, 
both at T = and T = 0.05, of the superfiuid density p s 
and the density p as functions of the chemical potential 
p for three different values of the onsite interaction Uq 
(= 4,5,6). [In both spinless and spin-1 cases, if U2 = 0, 
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TABLE I: The superfluid density p s , (F) 2 , and the compress- 
ibility k in the different phases of the spin-1 Bose- Hubbard 
model. The MI and NBL are really the same phase (see text). 
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FIG. 3: Plots of the variational free energy T (ground-state 
energy So for T = 0) as a function of ip for different values 
of p in the vicinity of the SF-MI transition for U% — and 
with = tpo = ipi = ip ( see text). The minima at ip = 
and ip 7^ correspond, respectively, to MI and SF phases, 
(a) The two minima at ip ^= merge into one minimum at 
ip = to yield the mean-field continuous SF-MI transition at 
T — as we increase p from 1.5 (bottom curve) to 2.3 (top 
curve) in steps of 0.2. Similar plots of T are given in (b) and 
(c). For < T, T develops three degenerate minima at the 
SF-MI boundary, indicating clearly the coexistence of SF and 
MI phases at a first-order boundary. This boundary can be 
crossed either (b) by changing p (from 0.7 to 1.5 in steps of 
0.2) at fixed T = 0.05 or (c) by changing T (0 to 0.1 in steps 
of 0.01) at fixed p = 2. 



the tip of the first lobe^iii lies at Uq c (p — 1) ~ 5.8 for 
T = 0.] Figure Q] (a) shows that Uq c (p) decreases as the 
T increases, i.e., the MI lobes grow at the expense of 
the SF phase. Figure [2] shows that p s goes to zero con- 
tinuously at the SF-MI transition, if T = 0, but with a 
jump if T = 0.05. Thus the SF-MI transition becomes a 
first-order transition at finite T and the zero-temperature 
SF-MI boundaries [Fig. Q](a)]are really lines of tricritical 
points; as the temperature is increased further, the first- 
order transition again becomes continuous at another tri- 
critical point. 

The first-order, SF-MI coexistence boundary is associ- 
ated with the three-degenerate-minima structure in the 
variational-free-energy plots shown in Fig. [3J To obtain 
these plots we use ip-i = ipo = ipi = ip that is one of 
the admissible solutions in the infinitely degenerate SF 
phase for the case J7 2 = 0; the infinite degeneracy of this 
phase, illustrated for T — in Fig. [1] (b), persists in our 
mean-field theory even if T > 0. Figure [3J shows plots 
of the variational free energy T (ground-state energy £q 
for T — 0) as a function of ip for different values of p, 
in the vicinity of the SF-MI transition; the minima at 
■0 = and ip ^ correspond, respectively, to MI and SF 
phases. At T = the SF-MI transition is continuous: 
this is reflected in the plots of £q in Fig. [3] (a), where, as 
we go from the SF to the MI phase by changing /i, two 
global minima with ip ^ merge to yield one minimum 
at ip = 0; precisely at the mean-field critical point we 
have a quartic minimum. For < T, J- develops three 
degenerate minima at the SF-MI boundary, indicating 
clearly the coexistence of SF and MI phases at a first- 
order boundary. This boundary can be crossed either by 
changing p at fixed T [Fig. [3] (b)] or by changing T at 
fixed p, [Fig. [3J (c)]. At sufficiently high temperatures 
this three-minima structure of T goes away at a tricrit- 
ical point at which the three minima coalesce to yield a 
sixth- order minimum. Beyond this tricritical point the 
SF-MI transition is continuous (second-order). 

Calculations such as those summarized in the plots of 
Fig. [3] help us to obtain the phase diagrams shown in 
Figs. H (a) - (d) for U = 12 and U 2 = 0. Let us begin 
with the p — T phase diagram shown in Fig. 0] (a). The 
MI phases [lobes in Fig. []] (a)] at T = evolve with- 
out any singularity into the normal Bose liquid (NBL) 
for T > 0. As we have emphasized earlier, MI and NBL 
phases are not distinct, but it is useful to think of a 
smooth crossover from one to the other; we define these 
crossover boundaries as the loci of points at which the 
compressibility k — k x = 0.02. The MI-NBL crossover 
boundaries (lines with filled triangles) are also shown in 
Figs. [4] (a) and (c). Islands of the SF phase appear in the 
fJi — T phase diagram; the first two of these are shown in 
Fig. 0] (a), where one is marked SF and the other, near 
the origin, is shown magnified in Figs. 0](b) and (c). The 
only difference between Figs. [4] (b) and (c) is that the 
latter shows the MI-NBL crossover boundary (line with 
triangles) and a line with stars, the locus of points in 
the SF phase at which the variational free energy T goes 
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FIG. 4: (a) Mean-field phase diagram in the fj, — T plane for Uo — 12 and U2 = 0. Lines with open (filled) circles represent 
first-order (continuous) SF-MI/NBL phase boundaries. First-order and continuous boundaries meet at tricritical points(TCP). 
The T = (T > 0) tricritical points are labeled TCP01, TCP02, etc. (TCP1, TCP2, etc. ). The line with triangles represents 
the crossover boundary between MI and NBL regions of the MI/NBL phase. The lower left corner of the phase diagram in (a) 
is enlarged in (b) and (c). The only difference between (b) and (c) is that the latter shows the MI-NBL crossover boundary 
(line with triangles) and a line with stars, the locus of points in the SF phase at which the variational free energy T goes from a 
curve with three minima to one with two minima; this line meets the SF-MI boundary at TCP1. (d) The density-temperature 
(p — T) version of part of the p — T phase diagram of (a) (without the MI-NBL crossover line); tie lines are used to hatch the 
two-phase regions in which SF and MI/NBL phases coexist (see text). 



from a curve with three minima to one with two min- 
ima. This line meets the SF-MI boundary at a tricriti- 
cal point labeled TCP1. Higher islands of the SF phase 
show analogous tricritical points labeled TCP2, TCP3, 
etc.; the SF-MI phase boundaries meet the T = axis at 
the zero-temperature tricritical points TCP01, TCP02, 
TCP03, etc. [Fig. gl(a)]. Figure|2(d) shows the density- 
temperature (j) — T) version of part of the /jl—T phase di- 
agram of Fig. [2(a) (without the MI-NBL crossover line); 
the first-order parts of the SF-MI boundaries now appear 
as regions of two-phase coexistence that are hatched with 
tie lines; the two-phase regions corresponding to the two 



lowermost SF-MI boundaries in Fig. [2 (a) are depicted; 
they end at the tricritical points TCP1 and TCP2 out of 
which emerge the continuous (second-order) SF-MI phase 
boundaries. We use the label MI/NBL since there is no 
strict distinction between MI and NBL phases for T > 0. 
Note that, in such a p — T phase diagram, the MI/NBL 
phases get pinched into exponentially small regions [e.g., 
in the vicinity of p = 1 in Fig. |2 (d)] as T — > and two 
zero-temperature tricrtical points get mapped onto each 
other [e.g., TCP01 and TCP02 in Fig. [2(d)]. 

We now investigate the case U2 ^ 0, so the onsite 
interaction between bosons depends on the spin. This 
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lifts some of the infinite degeneracy we encountered in 
the case U 2 = as can be seen directly at T = by 
comparing the pseudo-grayscale plots of £q in Figs. Q] 
(b), (c), and (d), for U 2 = 0, U 2 > 0, and U 2 < 0, 
respectively. 

If U 2 < there are four degenerate minima of, each 
corresponding to a ferromagnetic SF, with ip^i = ipi 
and tpo — iv^ii' The zero-temperature, mean-field, 
phase diagram for this case is shown in Fig. [51 it has the 
same topology as the phase diagram for the case U 2 = 
[Fig. [T] (a)]. We see that the MI phases have shrunk 
marginally and the SF-MI transitions are still continuous. 
The continuous nature of the T = 0, SF-MI transition is 
illustrated by the continuous variation of tp±ii' l Poi an d 
p s as functions of p in Fig. [51 The parameter (F) 2 , 
defined only in the SF phase, assumes the value 1, which 
confirms that we have a ferromagnetic SF phase in this 
case. The p — T phase diagram for the case U 2 < has 
the same topology as the U 2 = phase diagrams of Fig. 
UJ We do not show the p — T phase diagram for U 2 < 
since, for the parameters we use, namely, Uq = 12 and 
U 2 /Uq — —0.03, the phase boundaries are very close to 
those in Fig. [4j 
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FIG. 5: Mean-field phase diagram in the (p — Uo) plane for 
U2/U0 = —0.03 and T = 0. This has the same topology as the 
phase diagram for the case U2 = [Fig. [TJ (a) for T — 0] but 
the MI lobes have shrunk marginally; the SF-MI transitions 
are continuous. 

If U2 > there are four degenerate minima of the 
variational free energy T shown, e.g., at T — in the 
pseudo-grayscale plot of Fig. [T] Each one of these min- 
ima corresponds to a polar SF, with either = ipi 7^ 
and -00 = or vice versa as shown in the plots of %j>±\ and 



i/jo, versus p in Fig. [7] for U 2 /Uq = 0.03. The parameter 
(F), defined only in the SF phase, assumes the value 0, 
which also confirms that we have a polar SF phase. The 
zero-temperature, mean-field, phase diagram for this case 
is shown in Fig. [51 If the density p is equal to an odd 
integer (p = 1 is shown in Fig. [5J, this phase diagram 
has the same form as its counterpart for U 2 — [Fig. [T] 
(a)]. We see that the MI lobes expand marginally, and 
the SF-MI transitions are still continuous. However, if 
the density p is equal to an even number (p = 2 is shown 
in Fig. [SJ, the SF-MI transition becomes first-order and 
the MI phase is stable over a much wider region of param- 
eter space than for the case U 2 = 0: As we show in Fig. 
M for U 2 /U = 0.03, Uq — 7 and T — 0, p s and p vary 
continuously as functions of p at the SF-MI transition 
for p = 1 but discontinuously for p — 2; for comparison 
we also include the analogous plots for U 2 = 0. [We use 
Uq = 7 here, rather than Uo = 12, to compress the range 
of p, over which the SF-MI transitions occur.] 
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FIG. 6: Mean-field values of the superfluid order parameters 
$i=0_i and ipo and p a plotted as functions of p for Uo ~ 12, 
U2/U0 ~ —0.03, and T — 0; p s goes to zero continuously at 
the SF-MI transitions. The parameter {F) 2 , defined only in 
the SF phase, assumes the value 1, which confirms that we 
have a ferromagnetic SF phase in this case. 

For p = 2 in the MI phase, there are exactly two bosons 
localized per site and the total spin at every site can be 
either S — or S — 2. Since U 2 > there is an energy 
difference between the S = and 5 = 2 states, with a 
lower energy for the singlet state. To go from the MI to 
the SF phase, this singlet state has to be broken by sup- 
plying an energy ~ U 2 , which gives a rough estimate for 
the latent heat of this first-order transition if < T . This 
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requirement of a latent heat leads to the greater stability 
of the MI phases for even values of p relative to their 
counterparts for odd values of p. Thus the p = 2 MI lobe 
in Fig. [5] is substantially larger than the one for p = 1. 
The p — T phase diagram for the case U2 > 0, shown in 
Fig. HU](a), for U = 12 and U 2 /U = 0.03, has nearly the 
same form as the U2 — phase diagram of Fig. [J] The 
principal qualitative difference between these phase dia- 
grams is that, if t/2 > 0, there are no zero-temperature, 
tricritical points for the first-order boundaries associated 
with the MI lobes for even values of p; e.g., the tricritical 
point TCP03 in Fig. 0] has no counterpart in Fig. [TU] 

(a) . A quantitative comparison between these two phase 
diagrams is made in Fig. [TOJ(b); this shows that the two 
phase diagrams are nearly indistinguishable except for 
the first-order boundaries that link the zero-temperature, 
SF-MI transitions for even values of p with the tricriti- 
cal points directly above them (e.g., TCP3). In Fig. [TO] 

(b) the dashed line with open circles (open diamonds) is 
the first-order boundary for U2/U0 = 0.03 (U2 = 0); the 
region I between these lines lies in the MI (SF) phase 
if U2/U0 = 0.03 (U 2 = 0); the lines with filled triangles 
show the MI-NBL crossover as in Fig. @] (a) . Phase di- 
agrams such as Fig. [10] are obtained by calculating the 
order parameters ip a , and thence p and p Sl as functions 
of p at different temperatures. Representative plots are 
shown in Fig. [II] for U 2 /U = 0.3, U = 7, and T = 
and T = 0.05. 
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FIG. 8: Mean-field phase diagram in the (p — Uo) plane for 
U2/U0 = 0.03 and T = 0. The p = 1 MI lobe has the same 
form as its counterpart for U2 — [Fig. [1] (a)] . We see that 
the MI lobes expand marginally, and the SF-MI transitions 
are continuous (represented by a continuous line); but for 
p = 2 the SF-MI transition becomes first-order (represented 
by a dashed line) and the MI phase is stable over a much 
wider region of parameter space than for the case U2 = 0. 




FIG. 7: Mean-field values of the superfluid order parameters 
ipi — ip-i and ipo plotted as functions of p for Uo = 10, 
U2/U0 = 0.03, and T = 0. The SF phase has either V-i = 
ipi 7^ and ipo = or vice versa, which confirms that we have 
a polar SF in this case. 



IV. CONCLUSIONS 

We have carried out the most extensive study of 
the phase diagram of the spin-1 Bose Hubbard model 
so far by generalizing an intuitively appealing mean- 
field theory that has been used earlier for the spinlcss 
case. Our study yields both zero-temperature and finite- 
temperature phase diagrams for this model. Only T = 
phase diagrams had been obtained so far— 1 10 ' 11 ; so our 
elucidation of the finite-temperature properties of this 
model yields qualitatively new insights. We find, in 
particular, that several of the SF-MI transitions in this 
model are generically first order; at sufficiently high tem- 
peratures they become continuous via tricritical points. 
Tricritical points also abound at zero temperature since 
some, but not all, of the finite-temperature, first-order 
transitions become continuous as T — > 0. The resulting 
phase diagrams (Figs. [4land [T0]) are very rich and should 
provide a challenge for experimental studies, which we 
hope our work will stimulate. Experiments can study 
both the case U2 < 0, which can be realized possibly by 
using 87 Kb, and the case U% > 0, which can be realized by 
using 23 Na. Thus, in principle, both the phase diagrams 
of Figs. S] and [TO] could be obtained experimentally. Of 
course this will require good experimental control of both 
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the temperature and the density (or chemical potential) 
of the bosons. 

Our mean-field theory has been designed to investi- 
gate the relative stabilities of SF and MI/NBL phases. 
It has enough structure to unravel the differences be- 
tween polar and ferromagnetic supcfluids. However, our 
mean-field theory does not account for order parameters 
that can distinguish between different spin orderings in 
the MI phase, e.g., spin-singlet and nematic Mis, which 
have been investigated in the limit Uq — > oo by some 
group a 14 i 15 . The generalization of our mean-field theory 
to include such types of structures in the Ml phases of 
model (p} lies beyond the scope of this study but is an 
interesting challenge for further theoretical work. 



FIG. 9: Plots of p a and p as functions of p for T = 0, Uo = 7, 
and U2/U0 = 0.03 (filled circles) and U2 = (open circles). In 
the former case p 3 changes continuously at the SF-MI tran- 
sitions at the boundary of the p = 1 MI lobe (Fig. [5} but 
jumps at the first-order SF-MI transitions associated with the 
boundary of the p — 2 MI lobe. For U2 = only the p = 1 
lobe is encountered in this plot and the SF-MI transitions are 
continuous; p 3 shows a gentle minimum in the vicinity of the 
p = 2 MI lobe. 
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FIG. 10: (a) Mean-field phase diagram in the (p — T) plane for Uo = 12, U2/U0 ~ 0.03 showing first-order (open circles) 
and continuous (filled circles) transitions between SF and MI/NBL phases and tricritical points (TCP). This phase diagram is 
nearly the same as the U2 = phase diagram [Fig. 2], but there are no zero-temperature, tricritical points for the first-order 
boundaries associated with the MI lobes for even values of p: e.g., TCP03 in Fig. [4] has no counterpart here, (b) A quantitative 
comparison between these two phase diagrams shows that the two phase diagrams are nearly indistinguishable except for the 
first-order boundaries that link the zero-temperature, SF-MI transitions for even values of p with the tricritical points directly 
above them (TCP3 here); lines with open circles (open diamonds) denote the first-order boundaries for U2/U0 = 0.03 (U2 = 0); 
the region I between these lines lies in the MI (SF) phase if U2/U0 ~ 0.03 (U2 = 0); the lines with filled triangles show the 
MI-NBL crossover boundary. 
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FIG. 11: Representative plots of p 3 and p versus p for (a) T = 
and (b) T = 0.05 for U = 7, U2/U0 = 0.03 showing jumps 
at first-order SF-MI transitions; p a changes continuously at 
the T = 0, continuous SF-MI transition associated with the 
p = 1 MI lobe. 



